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Abstract 

The decays of a scalar particle, of either parity, into two photons or into two 
gravitons are evaluated. The effective interactions are of the form 4>FF, (pFF 
or (j)RR, (j)RR; in particular, Higgs boson decay into gravitons cannot be 
adduced to an interaction <j)R, as has been recently claimed. 

12.10.Gq, 12.10.Dm, 12.15.Ji, 12.25+e 



1 



Typeset using REVTp^X 



In a very recent paper, Srivastava and Widom jl| have claimed that the decay width of 
the Higgs meson into two gravitons is given by \/2G pm^ H / IQti . In their result, which they say 
stems from an effective interaction <j)R/ ((f)), Newton's constant disappears and gets replaced 
by the Fermi constant, leading to a large magnitude for the process. If their result were 
true, it would be counterintuitive to the notion that gravitational interactions are miniscule 
in particle physics and it would imply that the Higgs mesons disappears very quickly into 
a puff of gravitational plus electromagnetric radiation! In this paper we show that their 
result is not right and we derive the correct magnitudes for the decay amplitudes, be the 
decaying particle scalar or pseudoscalar; our faith in the weakness of induced gravitational 
(and electromagnetic) effects is happily restored. 

We consider parity-conserving decays of a + or 0~ particle into two photons (7) or 
into two gravitons (h). A simple helicity amplitude analysis shows that either process is 
governed by just one reduced helicity amplitude, (k,X;k',X'\S\0), where A = A' = 1 or 2, 
because of parity conservation. The same conclusion is reached by a covariant amplitude 
analysis, using the twin principles of gauge and general covariance. If e represents the 
external wavefunction of an incoming massless particle (7 or h), we may write down the 
unique couplings, 
Scalar — > j(k, X)j(k', A') 

Ccfryy = 9 e \(k)Ey{k )T pp {k } k ); T pp {k } k ) = k p k p — r/^pk ■ k , (1) 
Pseudoscalar — > "f(k, \)^{k' , A') 

Afe 77 = g 5 er(k)e*f(k')k a k'% pa p, (2) 

Scalar — > h(k, X)h(k', A') 

C*Kh = GeT^sT^T^k, k')T ua (k, k') + T up (k, k')T, a (k, k')}/2; (3) 
Pseudoscalar — > h(k, X)h{k' ', A') 

C^hh = G b sT{k)eT{k')[t ppaP k a k^T U(T {k, k') + (3 other perms)]/4; (4) 

Above, 7] stands for the Minkowski metric, about which we have expanded the gravitational 
metric field {g^ v = r]^ + nh^ v ) and with respect to which Lorentz scalar contractions are 
made. 

There are four significant points about eqs (1) - (4). Firstly, we note that insistence 
on current and stress tensor conservation means that each photon is accompanied by one 
power of momentum and each graviton comes with two powers of momentum, entirely in 
keeping with low-energy theorems ||. Secondly, we can encapsulate the gauge invariances 
by regarding the effective couplings above as corresponding to effective gauge interactions, 

g<j)F^F pu /4, gsfoe^F^F^/S, G<\>R pupa W vpG /2k 2 , G^e^BT a(i R paafi /An 2 , (5) 

where n 2 = 87rGNewton- Thirdly, the mass dimensions of g, g$ are M _1 and of G, G5 are M -3 . 
Fourthly and most importantly, the curvature tensor R has only two powers of derivatives, 
but of course contains all powers of the graviton field h, as befits a nonlinear theory. Thus it 
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is quite impossible for an interaction like <j)R to lead to two-graviton decay, as it is essential 
to conjure up four powers of momentum in the decay amplitude 0. In this respect we find 
the Srivastava-Widom claim extremely puzzling, notwithstanding that one is dealing with 
a spontaneously broken gauge symmetry and a vacuum expectation value which is closely 
tied to the Higgs boson field. 

It is straightforward to compute the decay rates ensuing from the interactions (1) to (4). 
Because gravitons and photons are massless the results are really quite simple in the end; if 
m denotes the mass of the decaying scalar, one finds the partial widths, 

r^ = \g\ 2 m 3 /QAn, = \g 5 \ 2 m 3 /QAn, r^ = 3|Gf m 7 /5127r, T^ hh = \G 5 \ 2 m 7 /512n. (6) 

The first and third of these are relevant to the Higgs boson. The only remaining question is 
the magnitude of the various coupling constants, g to G$. 

Because the decays are induced by quantum effects, they first arise at one-loop level. We 
may determine the couplings by applying the standard Feynman rules to the graphs depicted 
in Figure 1. In quoting the answers below, for a source field ip of mass circulating round 
the loop which is either scalar or spinor and of unit charge e, we should point out that certain 
graphs of Figure 1 are absent or give vanishing contributions; for example there is no Figure 
Id for photons, but it certainly arises for gravitons; similarly Figure lc does not occur for 
fermions undergoing electromagnetic decay but is present for gravitons and always exists 
for boson sources. As a purely technical point, we have applied dimensional regularization 
when calculating these diagrams, in order to respect gauge invariance. 

g — — ^ttt — t. / / dxdy 9(1 — x — y) (Axy — l)/(mi — m 2 xy); fermion loop (7a) 
2ir 2 Jo Jo 

£ 2 Qrh II 1"^ 

g = — ^p- / / dxdy 9(1 — x — y) xy/(m 2 h — m 2 xy); boson loop (7b) 
2ix 2 Jo Jo 
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J J dxdy 8(1 — x — y) / (m 2 ^ — m 2 xy); fermion loop only (8) 



G = / / dxdy 9(1 — x — y) xy(l — Axy) / (mL — m 2 xy)] fermion loop (9a) 

2ii l Jo Jo 

k 2 o^h ii 

G — / / dxdy 9(1 — x — y) x 2 y 2 /(ml — m 2 xy); boson loop (9b) 

2ii 2 Jo Jo 

G 5 = —LJbttt — t. I I dxdy 9(1 — x — y) xyl(mL — m 2 xy); fermion loop only. (10) 
2n z Jo Jo v 

Note that is dimensionless when ip is a spinor, but acquires dimensions of mass when 
ip is a scalar field. Anyhow, we see that all the integrals can be written as combinations of 
the basic integral, 
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In(v) = J q J dxdy6(l - x - y) {xy) n /(l - fixy); fi = m 2 /m^, (11) 

which can be expressed as sums of dilog functions and the like ||. As there is little to be 
gained by looking up tables of Spence functions, it turns out to be much easier to evaluate 
the I n numerically for input values of \i. 

Focussing on the Higgs boson now, and in respect of spin and spin 1/2 contributions 
to the decay widths r# 77 and r# ffff , we should remind ourselves that gu^ = -2my^ 
for scalar isosinglets and = —m^/v for quarks, where v = 246 GeV is the Higgs 

vacuum expectation value. Therefore the couplings entering the partial widths (6) equal g = 
—e 2 Ii(ii)/7r 2 v and G = —k 2 I 2 (ij)/ii 2 v from a boson loop as well as g = e 2 [I (fi)—4:Ii(fj J )}/n 2 v 
and G = — /t 2 [ii(/z) — 4/2 (/i)]/^ 2 ^ from a fermion loop; in particular, a massless loop particle 
(/i — > 00) reassuringly gives zero! Since the Higgs boson in the standard model couples 
just to the quarks/leptons and gauge bosons, apart from its self-interactions, we get the 
fermionic contributions, 



G- 



e 2 
2n 2 v 



-K 2 

2n 2 v 



E ^( J oK)-4J 1 (m 9 ))+^ Ul {m q )-4h{m g )) +Y,{Io(rn l )-4I 1 {m l )) 

q=u,c,t q=d,s,b l=e,^i,r 



E ^(AK)-4/ 2 K))+E i(/iK)-4i 2 K))+E(/iW-4/ 2 (m,)) 

—u,c,P q=d,s,b l=e,n,r 



(12) 



(13) 



To complete the calculation for electroweak theory one should add the contributions from 
the W and Z gauge bosons, including ghost field contributions as required by the gauge- 
fixing scheme in use. These have been computed previously |J, and are a factor of 21/4 
greater than the fermionic terms. But in any event we see that the gravitational decay rate 
is largely determined by the top quark (and the heavy gauge bosons), with a magnitude 
G ~ GNcwton^ 2 / 'vm 2 H that is fixed by the Newtonian constant. This enters eq. (6) and duly 
gives a miniscule result for the gravitational decay rate, of the order of 

?Hhh ~ m H (G Newton m 3 H /v) 2 /500 ~ m H x 1(T 75 . 

This is as it should be: the disappearance of the Higgs boson into gravitational radiation 
will never be observed. 
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FIGURE CAPTION 

Figure 1. One-loop contributions to Higgs decay into two massless gauge particles where 
the circulating loop particle can be scalar, spinor or vector. 
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FIGURES 




